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ABSTRACT 

Random  finite  set  provides  a  rigorous  foundation  for  opti¬ 
mal  Bayes  multi-target  filtering.  The  major  hurdle  faced  in 
Bayes  multi-target  filtering  is  the  inherent  computational  in¬ 
tractability.  Even  the  Probability  Hypothesis  Density  (PHD) 
filter,  which  propagates  only  the  first  moment  (or  PHD)  in¬ 
stead  of  the  full  multi-target  posterior,  still  involves  multiple 
integrals  with  no  closed  forms.  In  this  paper,  we  highlight 
the  relationship  between  Radon-Nikodym  derivative  and  set 
derivative  of  random  finite  sets  that  enables  a  Sequential 
Monte  Carlo  (SMC)  implementation  of  the  optimal  multi¬ 
target  filter.  In  addition,  a  generalised  SMC  method  to  im¬ 
plement  the  PHD  filter  is  also  presented.  The  SMC  PHD 
filter  has  an  attractive  feature-its  computational  complexity 
is  independent  of  the  (time-varying)  number  of  targets. 

Keywords'.  Multi-target  Tracking,  Optimal  Filtering,  Par¬ 
ticle  Methods,  Point  Processes,  Random  Sets,  Sequential 
Monte  Carlo. 

1.  INTRODUCTION 

Mahler’s  Finite  Set  Statistics  (FISST)  [8],  [9]  provides  a 
systematic  treatment  of  multi- sensor  multi-target  tracking 
using  random  set  theory.  In  FISST,  the  targets  are  treated 
as  a  single  meta  target,  called  the  multi-target  state,  and 
the  observations  are  treated  as  a  single  meta  measurement, 
called  the  multi-target  measurement.  The  multi-target  state 
and  multi-target  measurement  are  modeled  by  random  finite 
sets  whose  statistical  behaviour  are  characterised  by  belief 
functions.  The  key  to  rigorously  casting  the  multi-target  es¬ 
timation  problem  in  a  Bayesian  framework  is  the  notion  of 
belief  density,  based  on  set-derivatives  and  set-integrals. 

Analogous  to  the  single-target  case,  Bayes  multi-target 
tracking  propagates  the  multi-target  posterior  density  recur¬ 
sively  in  time  [9].  This  involves  the  evaluation  of  multi¬ 
ple  set-integrals  and  the  computational  intractability  is  far 
more  severe  than  its  single-target  counterpart.  The  Prob¬ 


ability  Hypothesis  Density  (PHD)  filter  [10],  which  only 
propagates  the  1  st  moment  (or  PHD)  of  the  multi-target  pos¬ 
terior,  requires  much  less  computational  power.  Unfortu¬ 
nately,  this  still  involves  multiple  integrals  that  in  gereral 
have  no  closed  form. 

Sequential  Monte  Carlo  (SMC)  methods  are  powerful 
tools  which  have  had  an  impact  on  optimal  (Bayesian)  fil¬ 
tering  [3].  The  key  idea  in  Monte  Carlo  methods  is  the  ap¬ 
proximation  of  integrals  by  random  samples  from  the  distri¬ 
bution  of  interest.  To  apply  SMC  methods  to  Bayes  multi¬ 
target  filtering  in  a  principled  way,  we  need  to  address  firstly, 
the  meaning  of  sampling  from  a  belief  density,  and  sec¬ 
ondly,  how  to  approximate  a  belief  density  from  these  sam¬ 
ples.  In  this  paper,  we  highlight  the  relationship  between 
Radon-Nikodym  derivative  and  set-derivative  of  random  fi¬ 
nite  sets  that  provides  answers  to  these  questions.  Con¬ 
sequently,  a  generic  SMC  implementation  of  the  optimal 
Bayes  multi-target  filter  can  be  proposed.  In  addition,  a 
novel  SMC  implementation  of  the  PHD  filter  which  offers 
a  much  cheaper  and  more  practical  alternative  is  also  pre¬ 
sented.  Both  algorithms  are  general  enough  to  capture  non¬ 
linear  non-Gaussian  dynamics. 

Section  2  of  the  paper  briefly  reviews  the  basics  of  ran¬ 
dom  finite  set,  optimal  Bayes  multi-target  tracking  and  de¬ 
scribes  the  PHD  or  1st  moment  filter.  Section  3  presents 
SMC  implementations  of  the  Bayes  multi-target  filter  and 
the  PHD  filter.  Simulation  results  are  presented  in  Section  4 
and  some  concluding  remarks  are  given  in  Section  5. 

2.  RANDOM  FINITE  SET  AND  BAYES 
MULTI-TARGET  FILTERING 

Multi-target  state  and  multi-target  measurement  at  time  k 
are  naturally  represented  as  finite  sets  Xk  and  Z^.  For  ex¬ 
ample,  if  at  time  k  there  are  M (k)  targets  Xk.i,  •  ■  Xk,M{k) 
in  the  state  space  Eg  then,  =  {xk.i,  ■  ■  ■  ,Xk^M{k)}  C 
Eg.  Similarly,  if  iV(fc)  observations  Zk.i,  ■  •  ■ ,  Zk^N(k)  in  the 
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observation  space  Eq  are  received  at  time  k,  then  = 
{zk,i,  ■  •  ■ ,  2fc,Ar(fc)}  C  Eo  where  some  of  the  N{k)  obser¬ 
vations  may  be  due  to  clutter.  Analogous  to  single-target 
system,  where  uncertainty  is  characterised  by  modelling  the 
state  and  measurement  by  random  vectors,  uncertainty  in 
a  multi-target  system  is  characterised  by  modelling  multi¬ 
target  state  and  multi-target  measurement  as  random  finite 
sets  (RFS)  Sfc  and  respectively.  A  formal  dehnition  of  a 
RFS  is  given  in  Section  2. 1 . 

Given  a  realisation  X^-i  of  Sfc_i,  the  multi-target  state 
at  time  k  can  be  modelled  by  the  RFS 

Ek  =  Sk{Xk-i)  U  Nk{Xk-i)  (1) 

where  Sk{Xk-i)  denotes  the  RFS  of  targets  that  have  sur¬ 
vived  at  time  k,  Nk{Xk-i)  =  BkiXk-i)  U  Ffc,  BkiXk-i) 
is  the  RFS  of  targets  spawned  from  X^-i  and  Ffc  is  the 
RFS  of  targets  that  appear  spontaneously  at  time  k.  The 
statistical  behaviour  of  the  RFS  5^  is  characterised  by  the 
conditional  probability  “density”  fk\k-i{Xk\Xk-i)  in  an 
analogous  fashion  to  the  Markov  transition  density  for  ran¬ 
dom  vector.  The  notion  of  probability  density  for  RFS  is 
formalised  in  Section  2.1. 

Similarly,  given  a  realisation  X^  of  5^,  the  multi-target 
measurement  can  be  modelled  by  the  RFS 

Sfc  =  Ok{Xk)  U  Ck{Xk)  (2) 

where  Qk{Xk)  denotes  the  RFS  of  measurements  generated 
by  Xk,  and  Ck{Xk)  denotes  the  RFS  of  clutter.  The  statisti¬ 
cal  behaviour  of  the  RFS  is  described  by  the  conditional 
probability  “density”  gk{Zk\Xk)  in  an  analogous  fashion  to 
the  likelihood  function  for  random  vector  observation. 

Let  Pfc  I  fc  ( Afc  I  Zo:fc )  denote  the  multi-target  posterior  “den¬ 
sity”.  Then,  the  optimal  multi-target  Bayes  filter  is  given  by 
the  recursion 

Pk\k-l{Xk\Zo:k-l) 

=  J  fk\k-iiXk\X)pk-iik-iiX\Zo..k-i)Ps{dX)(3) 

Pk\kiXk\Zo-.k) 

_  9k{Zk\Xk)pk\k-iiXk\Zo-.k-i) 

J  9k{Zk\X)pk\k-i{X\Zo.k-i)ps{dX) 

where  ps  is  a  dominating  measure  to  be  discussed  later  in 
Section  2. 1 .  The  main  difference  between  the  recursion  (3- 
4)  and  standard  clutter-free  single-target  filtering  is  that  Xk 
and  Zk  can  change  dimension  as  k  changes. 

2.1.  Random  Finite  Sets 

Given  a  closed  and  bounded  subset  E  of  R",  let  E{E)  de¬ 
note  the  collection  of  hnite  subsets  of  E.  A  random  hnite 
set  (RFS)  5  on  is  dehned  as  a  measurable  mapping 

S  :  G  ^  T{E). 


The  probability  measure  P  on  the  sample  space  G  induces 
a  probability  law  for  S,  which  can  be  specified  in  terms  of 
probability  distribution  or  belief  function.  The  probability 
distribution  is  defined  for  any  subset'  T  of  P{E)  by 

P^iT)  =  PE-\T)  =  Piliv  :  S(cu)  e  T}). 


The  probability  law  for  can  also  be  given  in  terms  of  the 
belief  function  [6],  [11],  dehned  for  any  subset  S  of  E 
by 

(3s{S)  =  P({cu  :  Eico)  C  S}). 

The  simplest  class  of  RFSs  are  the  Poisson  point  pro¬ 
cesses.  A  Poisson  point  process  is  a  RFS  T  with  the  prop¬ 
erty  that  for  any  k  disjoint  subsets  Si, ...,  Sk,  the  random 
variables  |T  n  S'! | , ...,  |T  n  S'fcj  are  independent  and  Pois¬ 
son.  Let  vr{S)  denote  the  mean  of  the  Poisson  random 
variable  |T  nS'|,  then  vy  dehnes  a  (dimensionless)  mea¬ 
sure  on  the  subsets  of  E,  and  is  called  the  intensity  measure 
of  T  [12].  The  probability  distribution  Py  is  given  by 


p-f(r)  = 

i=0 


v\{rnE^) 


where  denotes  the  ith  product  measure  of  vy  [4].  For 
any  subset  S  of  E,  let  Xk  {S)  denote  the  hyper- volume  of 
S  in  units  of  K.  The  density  of  vy  w.r.t.  Xk  is  called  the 
intensity  function  or  rate  of  T  and  has  units  of  K~^.  A 
Poisson  point  process  is  completely  characterised  by  its  in¬ 
tensity  measure  (or  equivalently  its  intensity  function). 

A  direct  extension  of  Bayesian  reasoning  to  multi-target 
systems  can  be  achieved  by  interpreting  the  density  ps  of  a 
RFS  5  as  the  Radon-Nikodym  derivative  of  the  correspond¬ 
ing  probability  distribution  P^  with  respect  to  an  appropri¬ 
ate  dominating  measure  p,  i.e 


PEiT)  =  j^pE{X)p{dX). 


The  dominating  measure  that  is  often  used  is  an  unnormalised 
distribution  of  a  Poisson  point  process  with  a  uniform  rate 
of  K~^  (intensity  measure  A  =  XkIK)  [4],  i.e. 


i=0 


(5) 


For  any  U  C  P{Es),  V  C  P{Eo)  let  Pk\k{id\Zo,k)  = 
P{Ek  G  U\Zij,k)  denote  the  posterior  probability  measure, 
Pk\k-i{U\Xk-i)  =  P{Ek  G  U\Xk-i)  and  Pk{V\Xk)  = 
P{X3k  G  V\Xk)  denote  the  conditional  probability  mea¬ 
sures  which  describes  the  multi-target  Markov  motion  and 
measurement  respectively.  Then  the  multi-target  posterior 
density  pk\k{-\ZQ,k),  transition  density  /fc|fc_i (-lAfe.i)  and 
likelihood  9k{-\Xk)  used  in  (3-4)  are  the  Radon-Nikodym 
derivatives  of  Pk\k{-\Zo:k),  Pk\k-i{-\Xk-i)  and  Pk{-\Xk) 
respectively. 

'where  appropriate,  subsets  are  assumed  to  be  Borel  measurable. 
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2.2.  Finite  Set  Statistics 

In  this  section,  key  concepts  in  finite  set  statistics  (FISST) 
and  its  relationship  with  conventional  point  process  theory 
are  highlighted.  In  particular,  we  outline  how  the  condi¬ 
tional  densities  /fe|fe_i(-|-)  and  gk{'\')  used  in  the  recursion 
(3-4)  can  be  systematically  constructed  from  the  underlying 
physical  model  of  the  sensors,  individual  target  dynamics, 
target  births  and  deaths  using  FISST. 

Individual  target  motion  in  a  multi-target  problem  is  of¬ 
ten  modelled  by  a  transition  density  on  the  state  space  Eg 
while  the  measurement  process  is  modelled  as  a  likelihood 
on  the  observation  space  Eg.  Consequently,  the  multi-target 
transition  density  and  likelihood  are  difficult  to  derive  since 
we  need  to  compute  Radon-Nikodym  derivatives  of  proba¬ 
bility  measures  on  the  subsets  of  E{Es)  and  E{Eo).  On  the 
other  hand,  belief  functions  are  defined  directly  on  the  sub¬ 
sets  of  Eg  and  Eo-  Hence,  models  for  multi-target  motion 
and  measurement  of  the  form 


Pk\k-i{S\Xk-i)  =  P{EkCS\Xk-i) 
f3k{T\Xk)  =  PiEkCT\Xk) 


can  be  systematically  constructed  (see  [9]).  However,  as  be¬ 
lief  functions  are  non-additive,  their  Radon-Nikodym  deriva¬ 
tives  (or  densities)  are  not  defined.  FISST  [6],  [9]  provides 
an  alternative  notion  of  density  for  a  RFS  based  on  the  set- 
derivative  of  its  belief  function  [6]. 

Let  B{E)  denote  the  collection  of  subsets  of  E.  The 
set-derivative  of  a  function  F  :  B{E)  [0,  oo)  at  a  point 
X  a  E\&&  mapping  {dF)x  :  B{E)  [0,  cxo)  defined  as  [6] 


idPUT) 


lim 


f(ruA^)-f(r) 

A/f  (Aa;) 


where  A^,  denotes  a  neighbourhood  of  x,  and  Ak(Ax)  is  its 
hyper-volume  in  units  of  K.  This  is  a  simplified  version  of 
the  complete  definition  given  in  [6] .  The  set  derivative  at  a 
finite  set  {xi , ...,  a:„}  is  defined  by  the  recursion 


where  {dF):^,  =  F  hy  convention. 

For  a  function  /  =  dF’(.)  (T),  the  set-integral  over  a  sub¬ 
set  S'  C  FI  is  defined  as  follows  [6],  [9] 


by  set-integrals  and  densities  are  replaced  by  set-derivatives 
of  corresponding  belief  functions  [9]. 

Note  that  {dF)x{9i)  has  unit  of  K~X\_  Hence,  the 
belief  density  evaluated  at  two  sets  Xi,  Xi  with  different 
number  of  elements  would  have  different  units.  This  lead 
to  problems  in  obtaining  the  maximum  a  posteriori  estimate 
since  the  units  of  the  posterior  density  are  not  commensu¬ 
rable  [9]. 

We  propose  a  unitless  set-derivative  by  replacing  A/c  in 
the  definition  of  set-derivatives  and  set-integrals  with  the  in¬ 
tensity  measure  A  of  a  unit  rate  Poisson  point  process.  This 
is  equivalent  to  introducing  a  unit  cancelling  constant.  It  is 
easy  to  see  that  the  generalised  fundamental  theorem  of  cal¬ 
culus  still  holds.  Moreover,  using  the  dominating  measure 
given  by  Eq.  (5),  we  have  for  any  U  = 

[  fiX)^^idX)=  [  fiX)SX.  (6) 

Ju  JS 

It  follows  from  Eq.  (6)  that 

(rf/3H)(.)(0)  =  dPs/dfi  =  pB.  (7) 

i.e.  the  set-derivative  of  the  belief  function  fds  is  the  den¬ 
sity  pb  of  the  corresponding  probability  measure  Pb  with 
respect  to  the  dominating  measure  p  (to  see  this,  note  that 
for  any  U  =  U“oS'%  fgPB(X)SX  =  f^pB(X)p(dX)  = 
Pb{U)  =  Ps{S)).  Consequently, 

fk\k-i{Xk\Xk-i)  =  {dPk\k-i{-\Xk-i))xki^)j 
9k{Zk\Xk)  =  {dPk\k{-\Xk))zk{^)- 

FISST  converts  the  construction  of  multi-target  densi¬ 
ties  from  multi-target  models  into  computing  set-derivatives 
of  belief  functions.  Procedures  for  analytically  differenti¬ 
ating  belief  functions  have  also  been  developed  [6],  [9]  to 
facilitate  the  task  for  practising  tracking  engineers. 

In  general,  the  multi- target  dynamic  model  (1-2)  yields 
the  following  multi-target  Markov  transition  and  likelihood 

fkik-iiXk\Xk-i) 

=  Y.  Sk\k-i(W\Xk-i)nk\k-i{Xk-W\Xk-iX^) 

WQXk 

gk{Zk\Xk)=  Y  Ok{W\Xk)ck{Zk-W\Xk).  (9) 

W<ZZk 


f{X)5X  =  Yl  [ .  f{{xi,--,Xi})Xxidxi...dxi). 

i=o  “'•S* 

Central  to  FISST  is  the  generalised  fundamental  theo¬ 
rem  of  calculus 

fix)  =  idF)xm  ^  FiS)  =  f  fiX)6X. 

Js 

In  the  FISST  framework,  the  optimal  Bayes  multi-target 
filter  has  the  same  form  as  (3-4)  where  integrals  are  replaced 


where  Sfc|fc-i  ('l-^fc-i)  is  the  density  of  the  RFS  SkiXk-i) 
of  surviving  targets,  nk\k-ii'\Xk-i)  is  the  density  of  the 
RFS  Nk{Xk-i)  of  new-born  targets,  9ki-\Xk)  is  the  density 
of  the  RFS  0fc(Afc)  of  target  generated  observations  and 
Cki'\Xk)  is  the  density  of  the  RFS  Ck{Xk)  of  false  alarms. 
Note  that  the  difference  operation  used  in  (8-9)  is  the  set  dif¬ 
ference.  The  reader  is  referred  to  [6],  [9]  for  details  on  how 
S/cIfc— 1  (' 1 ) j  '^k\k  —  li'\Xk—l)^  9ki'\Xk)  and  Cki'\Xk) 
can  be  derived  from  the  underlying  physical  model  of  the 
sensors,  individual  target  dynamics,  target  births  and  deaths. 
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2.3.  The  PHD  Filter 


3.  SEQUENTIAL  MONTE  CARLO 
IMPLEMENTATIONS 


The  probability  hypothesis  density  (PHD)  filter  is  a  1st  or¬ 
der  approximation  of  the  Bayes  multi-target  filter.  It  pro¬ 
vides  a  much  cheaper  alternative  by  propagating  the  1  st  mo¬ 
ment  in  stead  of  the  full  multi-target  posterior  [10]. 

A  finite  subset  X  €  X{E)  can  also  he  equivalently  rep¬ 
resented  hy  a  generalised  function  vvhere  5x  de¬ 

notes  the  Dirac  delta  function  centred  at  x.  Consequently, 
the  random  finite  set  S  can  also  be  represented  by  a  random 
density  ■  These  representations  are  commonly  used 

in  the  point  process  literature  [1],  [12]. 

Using  the  random  density  representation,  the  1st  mo¬ 
ment  or  (PHD)  D’e.  of  a  RFS  S  is  defined  hy 


De{x)  =  E 


(10) 


The  PHD  De  is  a  unique  function  (except  on  a  set  of  mea¬ 
sure  zero)  on  the  space  E.  Given  a  subset  S  C  E,  the  PHD 
measure  of  S,  i.e.  fg  DE(x)A(dx),  gives  the  expected  num¬ 
ber  of  elements  of  S  that  are  in  S.  The  peaks  of  the  PHD 
provide  estimates  for  the  elements  of  S. 

Let  7fc  denote  PHD  of  the  RFS  Ffc  of  targets  which 
appear  spontaneously;  6fc|fc-i  (■10  denote  the  PHD  of  the 
RFS  Hfc|fc_i({0)  spawned  by  a  target  with  previous  state 
0  efc|fc_i  (0  denote  the  probability  that  the  target  still  exist 
at  time  k  given  that  it  has  previous  state  fk\k-i  ( 'I  ')  de¬ 
note  the  transition  probability  density  of  individual  targets; 
gk{'  \  ■)  denote  the  likelihood  of  individual  targets;  Ck  de¬ 
note  clutter  probability  density;  Afc  denote  average  number 
of  Poisson  clutter  points  per  scan;  and  pu  denote  proba¬ 
bility  of  detection.  Define  the  PHD  prediction  and  update 
operators  $^1^-1,  respectively  as 


{^k\k-ia){x] 

1  =  j  <l>k\k-i{x,0oi{0Hd0  + 

lk{x), 

(11) 

(T'fcQ!)(x)  = 

v(x)  +  V 

a{x), 

(12) 

for  any  integrable  function  aon  Es,  where 

^k\k—l 

f k\k—l{^\  0  ^k\k—l{^ 

10, 

v{x) 

=  1-pd{x), 

^k,z{x) 

=  PD{x)gk{z\x), 

Kk{z) 

Let  Dk\k,  denote  the  PHD  of  the  multi-target  posterior 
Pk\k-  Assuming  that  the  RFS  involved  are  Poisson,  it  was 
shown  in  [10]  that  the  PHD  recursion  is  given  by 

Dk\k  =  {'^k  O  ^k\k-l)  {Dk-l\k-l)  ■  (13) 

Since  the  PHD  is  a  function  defined  on  the  space  where 
individual  targets  live,  its  propagation  requires  much  less 
computational  power  than  the  multi-target  posterior. 


The  propagation  (3-4)  of  the  multi-target  posterior  density 
recursively  in  time  involves  the  evaluation  of  multiple  set 
integrals  and  hence  the  computational  requirement  is  much 
more  intensive  than  single-target  filtering.  Sequential  Monte 
Carlo  (or  particle)  filtering  techniques  permits  recursive  prop¬ 
agation  of  the  set  of  weighted  random  samples  that  approx¬ 
imate  the  full  posterior  [3]. 

3.1.  SMC  implementation  of  the  optimal  Multi-target 
filter 

In  the  FIS  ST  framework,  what  does  it  mean  to  sample  from 
a  belief  density?  and  how  can  we  approximate  a  belief  den¬ 
sity  by  random  samples?.  The  answers  to  these  questions 
lie  in  Eq.  (7)  i.e.  the  belief  density  of  a  RFS  in  FISST  is  in¬ 
deed  a  density  of  the  corresponding  probability  distribution 
of  the  RFS  (see  Section  2.2).  Sampling  from  a  belief  den¬ 
sity  is  equivalent  to  sampling  from  the  probability  density 
and  approximating  a  belief  density  is  the  same  as  approxi¬ 
mating  a  probability  density. 

The  single-target  particle  filter  can  be  directly  gener¬ 
alised  to  the  multi-target  case.  In  the  multi-target  context, 
however,  each  particle  is  a  finite  set  and  the  particles  them¬ 
selves  can  thus  be  of  varying  dimensions. 

Assume  a  set  of  weighted  particles 
representing  the  multi-target  posterior  is  available. 

The  particle  filter  proceeds  to  approximate  the  multi-target 
posterior  at  time  fc  by  a  new  set  of  weighted  particles 

{w^k  \  as  follows 


Particle  Multi-target  Filter 

At  time  fc  >  1, 

Step  1:  Sampling  Step 

•  For  i  =  1,  ...,7V,  sample  (jfc  ^ 'I  and  set 


~(i)  (^^k  \ 

fk\k-l  { 

'4'’|4‘2i) 

•  Normalise  weights:  ZliLi  4*^  =  ^ 

Step  2:  Resampling  Step 

1, 

(  ( i\  ~7t\'I 

1  ^ 

•  Resample  7  4  j 

\  to  get  < 

1  i=l 

(4^’4’},^ 

The  importance  sampling  density  qk  {■\Xk-i,  is  a 
multi-target  density  and  x\^^  ~  9fe(  'I  Zk)  is  a  sam¬ 

ple  from  a  RFS  or  point  process.  Sampling  from  a  point 
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process  is  well  studied  see  for  example  [1],  [12]  and  refer¬ 
ences  therein. 

After  the  resampling  step,  an  optional  MCMC  step  can 
also  be  applied  to  increase  particle  diversity  [5].  Since  the 
particles  reside  in  spaces  of  different  dimensions,  a  reversible 
jump  MCMC  step  [7]  is  required.  Under  standard  assump¬ 
tions,  convergence  results  for  particle  filters  also  apply  to 
the  multi-target  case  [2]. 

From  the  multi-target  posterior  estimate  of  the 

target  set  at  time  k  can  be  given  by  the  expected  a  poste¬ 
riori  estimator  [6],  [9].  This  estimator  is  the  1st  moment 
(or  PHD)  Dk\k  (see  Section  2.3  Eq.  (10)).  Given  a  particle 

approximation  of  Pk\k^  the  particle  approx¬ 

imation  of  Dk\k  is  then 

Dk\k{x)  =  j 6y{x)Pk\k{.dX\ZQ,k) 
yex 
N 


and  compute  the  predicted  weights 

/_^S|  li)  — ^ ,  I  —  1,  -Lk-1 


Ui)  I  <Ik(^k  \^k-l’^k 


“fclfc-l 


(iFFT 

JkPk  ^  ^  1"^^} 


Step  2:  Update 

•  For  each  z  ^  ,  compute 

Jfc 

Dfc(U=  Y. 

j=l 

•  Fori  =  1,  ...,Lk-i  +  Jk,  update  weights 


~{i) 

K  = 


(4*’)+ 


Kk(z)  +  Ckiz) 


~(i) 

"fc|fc-r 


Step  3:  Resampling 

•  Compute  the  totai  mass  Nk\k  =  4^^’ 


to  get 


41 


{-(i)  , ^  ^k-l  +  Jk 

Rescaie  the  weights  Nk\k  to  get  jui^ ,  | 


k\k  j 

Lk 
i=l  ' 


The  main  practical  problem  with  the  multi-target  par¬ 
ticle  filter  is  the  need  to  perform  importance  sampling  in 
very  high  dimensional  spaces  if  many  targets  are  present. 
Moreover,  it  can  be  difficult  to  find  an  efficient  importance 
density  and  the  choice  of  a  naive  importance  density  like 
qk{-\xj;^l^,Zk)  =  /fc|fc-i(-|  will  typically  lead  to 

an  algorithm  whose  efficiency  decreases  exponentially  with 
the  number  of  targets  for  a  fixed  number  of  particles. 

3.2.  SMC  implementation  of  the  PHD  filter 

The  PHD  filter  is  a  cheaper  alternative.  However,  direct  ap¬ 
plication  of  standard  SMC  methods  to  propagate  the  multi¬ 
target  PHD  would  fail  because  firstly,  the  PHD  is  not  a  prob¬ 
ability  density  function;  and  secondly,  it  is  not  a  standard 
Bayes  recursion.  In  this  section  we  summarise  the  SMC 
implementation  of  the  PHD  filter  proposed  in  [13]. 

For  any  fc  >  0,  let  ak  =  {Wk  \xk  ^}i=i  denote  a  parti¬ 
cle  approximation  of  Dk\k-  Using  the  PHD  recursion,  a  par¬ 
ticle  approximation  of  the  PHD  at  time  step  k  >  0  can  be 
obtained  from  a  particle  approximation  at  the  previous  time 
step  by  the  following  procedure  (see  [13]  for  the  derivation) 


SMC  PHD  filter 

At  time  A;  >  1, 

Step  1:  Prediction 

•  Sample 

~(i)  j  )  *  =  !>  ••■)  i-fc-l 

\  Pk('\ Zk)  X  =  Lk-l  +  ■■■,Lk-l  +  Jk 


Note  that  the  computational  complexity  of  this  algo¬ 
rithm  is  independent  of  the  (time-varying)  number  of  tar¬ 
gets  if  we  fix  the  number  of  particles  Lk  =  L  for  all  k. 
However,  this  can  result  in  instances  where  the  number  of 
particles  are  not  sufficient  to  resolve  a  large  number  of  tar¬ 
gets  while  at  other  times  we  may  have  an  excess  of  parti¬ 
cles  for  a  small  number  of  targets.  On  the  other  hand,  if 
Lk  =  Lk-1  +  Jk,  then  Lk  increases  over  time  even  if  the 
number  of  targets  does  not.  This  is  very  inefficient,  since 
computational  resource  is  wasted  in  exploring  regions  of  the 
state  space  where  there  are  no  targets.  It  would  be  computa¬ 
tionally  more  efficient  to  adaptively  allocate  approximately 
say  p  particles  per  target  at  each  time  step. 

Since  the  expected  number  of  targets  J  Dk\k{^\Z[).k)d^ 
can  be  estimated  by  Nk\k  —  ^4’  h  is  natu¬ 

ral  to  have  the  number  of  particles  Lk  =  pNk\k-  This 
can  be  achieved  by  resampling  Lk  =  pNk\k  particles  from 

redistributing  the  total  mass  Nk\k 
among  the  Lk  resampled  particles. 

In  the  resampling  step  of  the  SMC  PHD  filter,  the  new 
weights  normalised  to  1  (as  in  conven¬ 

tional  particle  filter)  but  sum  to  Nk\k  instead.  Similarly  to 
the  standard  case,  each  particle  x}.  is  copied  times  un¬ 
der  the  constraint  ~  Tfc  to  obtain  {x[,*4i4i' 

The  (random)  resampling  mechanism  is  chosen  such  that 
E;[iv4]  =  Lka^k^  where  >  0,  “4  =  1  is 

a  sequence  of  weights  set  by  the  user.  However,  the  new 
weights  are  set  to  oc  Wk'^ /a^k  "'i’-i^  X4i  ^  =  -^fc|fc 
instead  of  XXi  ^  =  1.  Typically,  a^k  =  Wk'’ /Nk\k- 
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4.  SIMULATIONS 

For  visualisation  purposes,  consider  a  one-dimensional  sce¬ 
nario  with  an  unknown  and  time  varying  number  of  targets 
observed  in  clutter.  The  states  of  the  targets  consist  of  1- 
D  position  and  1-D  velocity,  while  only  position  measure¬ 
ments  are  available.  Without  loss  of  generality,  we  con¬ 
sider  targets  with  linear  Gaussian  dynamics.  We  assume  a 
Poisson  model  for  spontaneous  target  birth  with  intensity 
0.2N  ( ‘I  0, 1).  Each  existing  a  target  has  a  (state  indepen¬ 
dent)  probability  of  survival  e  =  0.8.  The  clutter  process  is 
Poisson  with  uniform  intensity  over  the  region  [—100;  100] 
and  has  an  average  rate  of  10.  Figure  1  shows  the  tracks  with 
clutter  on  the  position  measurements  and  Figure  2  plots  the 
PHD  of  position  against  time.  Observe  from  Figure  2  that 
the  PHD  filter  shows  surprisingly  good  performance,  even 
very  short  tracks  are  picked  up  among  clutter. 


Fig.  1.  birth  rate  0.2,  death  rate  0.2,  clutter  rate  10 


Fig.  2.  PHD  of  position 


5.  CONCLUSION 

In  this  paper,  we  have  highlighted  the  relationship  between 
Radon-Nikodym  derivative  and  set  derivative  of  random  fi¬ 


nite  sets.  This  allows  the  derivation  of  a  generic  particle  fil¬ 
ter  for  estimating  an  evolving  random  finite  set.  In  addition, 
a  sequential  Monte  Carlo  implementation  of  the  probabil¬ 
ity  hypothesis  density  filter  for  multi-target  tracking  is  also 
given.  With  a  computational  complexity  that  is  independent 
of  the  (time-varying)  number  of  targets,  this  offers  a  very 
practical  alternative  to  the  particle  random  set  filter. 
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